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ABSTRACT 

Context. Short timescale flux density variations in flat spectrum radio sources are often explained by the scattering of radio waves 
in the turbulent ionized interstellar medium of the Milky Way. One of the most convincing observational arguments in favor of 
this is the annual modulation of the variability timescale caused by the orbital motion of Earth around the Sun. Jl 128+5925 is a 
recently discovered IDV source with a possible annual modulation in its variability timescale. New observations suggest a change in 
its variability characteristics. 

Aims. We search for changes either in source structure or in the interstellar medium that can explain the variations in the IDV pattern 
of Jl 128+5925. Using Very Long Baseline Interferometry (VLBI), we study a possible relation between source orientation on the 
sky and anisotropy angle seen in the annual modulation. Additionally, obtaining a source size estimate from VLBI data enables us to 
calculate the distance to the screen causing the variations in Jl 128+5925. 

Methods. We observed the source in six consecutive epochs with the Very Long Baseline Array (VLBA) at three frequencies, 5 GHz, 
8 GHz, and 15 GHz in total intensity and polarization. This data are combined with our densely time-sampled flux-density monitoring 
performed with the radio telescopes at Effelsberg (Germany) and at Urumqi (China). 

Results. The VLBA observations detected an east-west oriented core-jet structure with no significant motion in its jet. The expansion 
of the VLBI core leads to an estimate of mild relativistic speed (2.5 c ± 1.4 c). The position angle of the VLBI jet agrees with the 
angle of anisotropy derived from the annual modulation model. No significant long-term structural changes were observed with VLBI 
on mas-scales, although, the VLBI core-size expansion offers a possible explanation of the observed decrease in the strength of IDV. 
VLBI polarimetry measured significant changes in the electric vector position angle (EVPA) and rotation measure (RM) of the core 
and jet. Part of the observed RM variability could be attributed to a scattering screen (37 pc distance) that covers the source (core and 
jet) and may be responsible for the IDV. Superposition of polarized subcomponents below the angular resolution limit may also affect 
the observed RM. 
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1. Introduction tic variability time-scale is inversely proportional to the rela- 

. . . . . . tive velocity between the observer and the scattering medium. 

Intaday variability (IDV, |Heeschen et al.| U283 IWitzel et al.| The so -called "annual modulation" of the variability timescale 

1 19861) is the rapid variation observed usually at cm wave- reflects the systematic variation in the relative velocity vec- 



dKellermann & Pauliny-Toth|[l969l). In this source-intrinsic in- 2001: lJauncev & Macquart 2001), PKS 15 19 -273 (iJauncev et al. 
terpretation, relativistic boosting with large Doppler factors of 120031 ICarter et all 120091). PKS 1257- 326 dBignaU et all 12006). 
» 50 would be required to reduce these brightness temperatures an( j PKS1622-253 ( Carter et al 112009b 



» 50 would be required to reduce these brightness temperatures an( j PKS1622-253 (Carter et al 2009) 
to the inverse Compton limit of 10 12 K. However, such high val- 
ues of Doppler factors were not observed in kinematic studies of Unfortunately, the detection of an annual cycle in IDV 
quasars and blazars. sources is not always straightforward: in some IDV sources, 
One solution to this controversy is the source-extrinsic the- a previous claim of an annual v ariability cycle cou l d not 
ory of IDV. This explains that the rapid variations are caused dur- be confirmed in subsequent year dKedziora-Chudczerl [2006, 
ing the propagation of radio waves. The variations are caused by PKS 0405-385), and in one source at least ceased completely 
interstellar scintillation (ISS) in the intervening t urbulent, ion- dFuhrmann et al.ll2.Q02i B0917+624). Since the ISS phenomena 
ized i nterstellar medium of the Milky Way (e.g., iRickett et al.l requires a source size of the order or smaller than the scattering 
2006, and references therein). In this scenario, the characteris- size of the screen, possible source intrinsic size variations (e.g., 
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expansion or blending effects of the scintillating component(s)) 
could explain these variations in the I DV pattern or even its 
disappearance (e.g., iRraus et ail [1999: Krichbaum e t alj 120021; 
Bern hart et al.ll200q) . Alternatively, changes in the IDV charac- 
teristi cs may also be caused by changes in the scattering plasma 
(e.g., lKedziora-Chudczerl l2006). One way of distinguishing be- 
tween these different possibilities, is to perform a polarimetric 
VLBI monitoring at different frequencies. 

Jl 128+5925 (hereafter Jl 128+592) is an intraday variable 
source. Densely time-sampled flux density monitoring observa- 
tions of the source have been carried out with the Effelsberg 100- 
meter (MPIfR, Germany) and the Urumqi 25 -meter (China) ra- 
dio telescopes at 5 GHz over the past four years. In most of these 
observations, Jl 128+592 showed pronounced variability with 
peak-to-trough amplitudes exceeding 20%. The measured char- 
acteristic variability timescales were differ significantly from 
epoch to epoch, ranging from 0.2 to 1.6 days. The changes in 
the variability timescales are indicative of an an nual cycl e. 

After the discovery of IDV in Jl 128+592. IWu~eTaT] d2008l) 
performed optical observation of the source. Jl 128+592 did not 
show r apid variability i n the optical bands. Additional obser- 
vations dWu et all 12009) confirmed the optical quietness of the 
source. The different variability behavior in optical and radio 
regimes is expected if the IDV in radio wavelengths caused by 
interstellar scintillation. 

In the case of J 1128+592, we observed in late 2007 and 
early 2008 a strong decrease in the IDV amplitudes. To find 
out whether source-intrinsic changes lead to a change in the 
variability behavior, we performed 6 epochs of Very Long 
Baseline Array (VLBA) observations between 2007 July and 
2008 February. In this paper, we report on an analysis of these 
VLBA observations. 

The paper is organized as follows. In Sect. [2] we provide a 
brief summary of the single dish observations and the annual 
modulation model. In Sect. [3] we describe the VLBA observa- 
tions and the data reduction procedure. In Sect. |H we present 
the total intensity results of the VLBA observations, and in Sect. 
14.11 we discuss the IDV properties of J 1 128+592 in the light of 
the VLBA observations, and in Sect. [5] we present the polarized 
intensity results of the observations, which we discuss in Sect. 
15.11 Finally, in Sect.[6]we summarize our findings. 

The following cosmological parameters were used 
throughout this paper: Hq = 71 kms^'Mpc -1 , Qmatter = 
0.27, and Qyac = 0.73. At the the redshift (z = 1.795) of 
the source, 1 mas corresponds to a spatial scale of 8.5 pc. An 
apparent angular separation rate of 0. 1 mas/yr corresponds to an 
apparent velocity of 7.8 c. 



2. Summary of the single-dish observations and the 
annual modulation model 

Since 2004 December, Jl 128+592 has been observed 28 times 
at 5 GHz with either the Effelsberg 100-meter radio telescope 
(in 10 epochs) or the Urum qi 25-meter radio te le scope (in 18 
epochs). In previous papers (Gaba nyi et al.l2 007a b. 2008), parts 
of these observations were analyzed and discussed. Regarding 
the observational techniques, data reduction, and timescale de- 
termination, we refer to those papers and references therein. 

Two 5 GHz light curves of Jl 128+592 are shown in Fig. 
Q] The variability is pronounced, the peak-to-trough variability 
amplitude reaching 25 %. The two light curves exhibit different 
variability timescales of ~ 0.2 day in February 2006 and more 
than 1.5 day in April 2006. We proposed that the changes in the 
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Fig. 1. 5 GHz light curves of Jl 128+592, observ ed in February 
2006 ( top panel) and April 2006 (bottom panel) dGabanvi et alj 
l2007bl) . Observations are performed with the Effelsberg antenna 
(gray symbols) and the Urumqi antenna (black symbols). In the 
lower parts of both panels, the light curve of a secondary cali- 
brator source (B0836+710) is displayed for comparison. 



characteristic variability timescale are caused by annual modu- 
lation. In this case, the variability timescale is proportional to the 
scattering length-scale, and inversely proportional to the relative 
velocity between the scattering screen and the observer. 

However, the simple isotropic annual modulation model 
could not reproduce all observ ed changes in the variability 
timescale (Ga banyi et alj|2007b|) . In a more detailed scenario, 
the scintillation patter n could be more acc urately described by 
an anisotropic model (Bign aTTet al.l l2006b . Here, the variabil- 
ity timescale also depends on both the ellipticity of the scin- 
tillation pattern and the direction in which the relative veloc- 
ity vector "cuts through" the elliptical scintillation pattern (and 
not just the absolute value of the velocity vector). We note that 
this model does not infer the origin of the anisotropy, be it 
source intrinsic or induced by the scattering medium. A simi- 
lar anisotropic model was used to explain the sea sonal cycles 
in J1 819+3845 (iDennett-Thorpe & de Bruvr3 l2003h. PKS1257- 
326 (Big nalletalJ 120061). an d in both PKSB1519-273 and 
PKSB1622-253 dCartSLStaD EQOl). We therefore fit t ed the 
anisotropic annual modulation model of Big nall et al.l {2006) 
to our measurements with the following model parameters: the 
components of the screen velocity vector (in right ascension 
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Fig. 2. The variability timescales of the IDV of J 1 128+592, mea- 
sured at different epochs. The dashed black line shows the fitted 
annual modulation model. Circles represent the data points used 
in the fitting. Open circles represent data from the Urumqi tele- 
scope, filled circles stand for Effelsberg measurements. Vertical 
lines mark the epochs of the VLBA observations. 
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Fig. 3. Flux density of Jl 128+592 as measured with single dish 
telescopes (open circles) and by the VLBA (filled gray dia- 
monds) at 5 GHz. Black dots represent the standard deviation of 
the flux density variations of Jl 128+592. The standard deviation 
values are multiplied by a factor of 10 for a clearer display. 



Table 1. Details of the VLBA observations of Jl 128+592. 



v a , and declination direction Vg), the scattering length scale (s, 
which depends on the screen distance and the scattering size), 
the position angle of anisotropy (/?), and the axial ratio of the 
anisotropy (r). 

As the monitoring observations continued and we gathered 
more data, we were able to refine the annual modulation model. 
For all available data acquired until 2008 April, we obtain the 
following parameters from the fits: v a = (3 + 4)km/s, v$ 
( 1 1 + 2)km/s, s = (1.3 + 0.3) • 10 6 km, r = 3.4 + 0.8, and 
/3 = -98° ± 5°. The angle is measured north through east. The 
measured timescales and the fitted annual modulation model are 
displayed in Fig. [2] 

During the 4 years of monitoring, the variability amplitude 
decreased. In Fig. [3] we plot the total flux density at 5 GHz 
versus time. We also plot the standard deviation of the individ- 
ual flux density measurements in a single observing campaign, 
which is a measure of the "strength of the IDV". A value cor- 
responding to ten times the standard deviation is displayed as 
filled circles in Fig. [3] We performed 6 epochs of VLBA obser- 
vations during 2007 and 2008, to investigate the physical reasons 
behind the decreasing variability strength and search for related 
structural variations in the source. 

Comparing the average flux densities of these observations, 
we found that the source experienced significant long-term flux 
density change (see Fig. [3). The flux density of Jl 128+592 in- 
creased by ~ 27 % until the end of 2006, afterward the flux den- 
sity continuously decreased until October 2007 by ~ 63%. In 
parallel to the decreasing total flux density, the strength of the 
variability (standard deviation given by filled circles in Fig. 
decreased in a similar way. We note that during the phase of 
brightening (before 2006), Jl 128+592 was not monitored in a 
systematic way. Therefore, owing to differences in the time sam- 
pling and duration of the individual measurements, the error bars 
and the scatter in the measured standard deviation are larger for 
the first couple of measurements (for a discussion of the long- 
term flux density variability, see Sect.|4]l. 
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2004.047 
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all 




2004.211 
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all 
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all 




2004.545 
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all 




2004.706 
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all except KP 




2004.884 
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all 
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2007.540 


5/8/15 


all, except HN at 15 GHz 
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2007.661 


5/8/15 


all, except HN 


C 


2007.773 


5/8/15 


all, except SC 


D 


2007.907 


5/8/15 


all, except SC and PT 


E 


2008.025 


5/8/15 


all 


F 


2008.137 


5/8/15 


all 



3. VLBA observations and data reduction 

Our first VLBA observation took place in July 2007, and 5 sub- 
sequent epochs were followed by a 4 to 6 week separation. The 
last VLBA observation was performed in February 2008. Each 
observation had a duration of 6 hrs and was performed at 5 GHz, 
8 GHz, and 15 GHz with dual polarization. The VLBI data were 
recorded with a data-rate of 256 Mbit/s. The sources 4C39.25 
and J 103 8 +05 12 were included as calibrators. In addition to 
the data from our 6 VLBA experiments, we also searched the 
VLBA archive for previous observations of Jl 128+592. We (re)- 
analyzed six epochs of 5 GHz observations closest int time to 
our monitoring program. J 1 128+592 was used as a calibrator in 
these archival observations. We unfortunately did not find 5 GHz 
data overlapping with our single-dish campaign. The details of 
all these observations are summarized in TableQ] In the table, the 
first 6 epochs are 5 GHz observations from the VLBA archive. 
The last 6 epochs are our multi-frequency observations. 

After correlation at the VLBA correlator in Socorro 
(NRAO), the data were calibrated and fringe-fitted using the 
standard procedures within the AIPS (Astronomical Image 
Processing System) software package. The post-processing of 
the data included the usual steps of editing, phase- and amplitude 
self-calibration, and imaging. These tasks were performed with 
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the Caltech DIFMAP program package. The Caltech DIFMAP 
package was also used to fit the uv-data. We fitted the source 
structure in every epoch with circular Gaussian components. We 
tested fitting with elliptical Gaussians, however we did not ob- 
tain stable and physically meaningful results for every epoch 
(e.g., elliptical Gaussian components became degenerate). For 
consistency and to follow the same reduction and analyzing steps 
in every epoch, we parametrized the brightness distribution of 
J 1128+592 by fitting only circular Gaussian components to the 
visibilities. 

Calibration of the polarized data involved the correction for 
parallactic angle, determination of both right-left multiband and 
single-band delay corrections, feed D-term calibration, and the 
absolute calibration of the electric vector position angle (EVPA). 
For D-term calibration, J1038+0512 was used. Its total intensity 
and polarization structure is simple and point-like, compared to 
our other calibrator, 4C39.25. The EVPA calibration was com- 
pleted by comparing the averaged EVPA of 4C39.25 measured 
in the VLBA images and those values given by the Very Long 
Array (VLA) polarization monitoring prograrrQ at 5 GHz and 
8 GHz. To calibrate EVPA at 15 GHz, we used the 14.5 GHz 
measurements of the Universi ty of Michigan R adio Astronomy 
Observatory (UMRAO, e.g.. lAller et all 120031) . The UMRAO 
observations are performed at three frequencies, 5 GHz, 8 GHz, 
and 14.5 GHz. The VLA and UMRAO results are all consis- 
tent, apart from one corresponding to the first epoch, for which 
a difference of ~ 100° between the 8 GHz EVPA of 4C39.25 
was measured. The observations at the different sites were per- 
formed within two days. The measured polarized flux density 
differed significantly, the VLA flux density being twice as high 
as the UMRAO one. Since the amount of polarized flux den- 
sity measured in our VLBA maps agreed with that measured for 
the UMRAO data, we decided to use the EVPA reported by the 
UMRAO. 

4. Results of the VLBA observations - total intensity 

The aperture synthesis imaging of a non-stationary source, 
which varies during the time of the observation, leads to a degra- 
dation in the reconstructed image. This amount of degradation 
depends on the variability amplitude and the timescale of the 
variability with respect to the total observing time. Simulations 
showed that a smooth change in 80 % in the source flux density 
over an observing time of 12hrs leads to residual side-lobes and 
a limitation of the d ynamic range i n the VLBI map of the order 
of a few hundred dHummellll987h . In the case of Jl 128+592, 
the maximum peak-to-trough variability amplitude during our 6 
VLBA experiments was ~ 16% (Dec. 2007), but is typically 
lower and of the order of ~ 10%. The variability timescales 
were in the range 0.3 - 1.3 days, which is extremely long com- 
pared to the 6hr duration of our VLBA experiments. We ex- 
amined our VLBI images and we found no strong signs of im- 
age degradation (excessive rms, side-lobes, or symmetric struc- 
tures indicative of calibration errors). At 5 GHz, where the IDV 
is strongest, we subdivided a VLBA observation into two 3 hr 
segments and imaged them separately (IDV decreases towards 
higher frequencies, resulting in much lower image degradation 
at 15 GHz). The resulting small differences in the flux densities 
of the VLBI components (< lOmJy) are within the errors intro- 
duced by the a-priori flux density calibration, which is accurate 
at the 5-10% level, and the accuracy achieved by Gaussian 

1 http : //www . via . nrao . edu/astro/calib/polar 



Table 2. Flux densities of the fitted Gaussian model components 
in the six archival VLBA epochs. 



Epoch 


Icore 


[mJy] 


/ jet [mJy] 


2004.047 


374 


±25 


125 ± 12 


2004.211 


378 


±40 


23 ± 2 


2004.386 


398 


±80 


116 ±80 


2004.545 


297 


±30 


144 ± 14 


2004.706 


546 


±35 


57 ±33 


2004.884 


434 


±20 


1 12 ± 10 



model fitting, which is mainly limited by the uv-coverage. We 
therefore conclude that the IDV, which was present during our 
VLBI observations, did not affect the resultant images signifi- 
cantly, nor could it cause largely wrong estimates of the VLBI 
component flux densities derived from Gaussian model fits. 

The VLBA maps (Fig. [4] Fig. [5] and Fig.[6]i show the east- 
west oriented two-component structure of Jl 128+592, which 
is clearly resolved at 8 GHz and 15 GHz. At 5 GHz, the to- 
tal intensity maps show a point-like source at first glance, but 
after model-fitting with circular Gaussian components, a two- 
component model is found to provide a significantly better fit to 
the data in every epoch. 

The brighter and more compact component has a flat spec- 
tral index between the 5 GHz and 15 GHfl Its average spectral 
index during the six epochs is -0.06 ± 0.04. The western com- 
ponent has a steep spectrum with an average spectral index of 
-1.0 + 0.1. Thus, we can conclude that the brighter component 
is most probably the compact core, while the westward feature 
is an optically thin jet component. 

The flux densities of the jet and core components at 5 GHz, 
8 GHz, and 15 GHz are listed in Tables [3] El and|5j respectively. 
The summed flux density of the two fitted components at 5 GHz 
agrees with the mean flux densities measured by the Urumqi 
telescope during that period (Fig.[3]and[7]l. More than 90% of 
the flux density measured by a single dish is recovered by our 
VLBA observations at 5 GHz. 

The archival 5 GHz VLBA data also shows an east-west ori- 
ented core-jet structure. However, the summed flux density of 
the model components are larger in every epoch than in our 
VLBA observations. This is in agreement with the long-term flux 
density changes revealed by the single-dish observations (Fig. 

According to the archival VLBA data (Table|2|, the flux den- 
sity of the jet feature is very similar to the results of our VLBA 
observations. But the core component is brighter than in our ob- 
servations in 2007 and 2008, which suggests that the core was 
responsible for the brightening observed in the single dish mon- 
itoring. 

According to our VLBA observations, the size of the fitted 
circular core component is 6$ qj-[ z = 0.2 + 0.05 mas at 5 GHz, 
6 8 GHz = (0.08 ± 0.02) mas at 8 GHz, and 6> 15 GHz ~ 0.15 mas 
at 15 GHz. In the standard jet model (e.g. Bfandford & Konigl 
1 1979b . the base of the jet is assumed to have a constant bright- 
ness temperature of 7g ~ S /(v 2 8 2 ) ~ constant (where S is the 
flux density, v the frequency, and 9 is the FWHM source size). 
Thus, the size of a flat spectrum source can be expected to be 
inversely proportional to the frequency. However, the size of the 
flat spectrum core of Jl 128+592 does not follow this relation. 
The 15 GHz size appears excessively large. This might indicate 
that the core component is a blend of the core and a new feature 



2 The spectral index is defined as S ~ V", where v is the frequency 
and 5 is the flux density. 
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Fig. 4. 5 GHz VLBA maps of Jl 128+592 from the six observing epochs. Contours show total intensity, lines show the strength and 
direction of the polarized intensity. Contours are in percent of the peak flux and increase by factors of two. The peak flux intensities 
are 312, 299, 295, 292, 277, and 300mJy/beam. The beam-size is 2 mas x 1 mas and is shown at the bottom left corner of each 
image. The 1 mas length of the superimposed polarization vectors corresponds to 1 mJy/beam. 



that is not (yet) resolved at 15 GHz at the time of our observa- 
tions. 

From the six reanalyzed archival epochs, we obtain an av- 
erage core size of ~ (0.07 + 0.02) mas. This is less than half 
of the value that we measured in our VLBA observations. 
Unfortunately, we do not have VLBI data acquired simultane- 
ously to the time when the source was in its most variable IDV 
state. The available data, however, indicate that the core might 
have been more compact at that time. Since a larger size of the 
scintillating compo nent(s) quenches the variations more (e.g., 
Rick ett et alj 1 19951 and references therein), this increase may 
naturally explain the observed decrease in the variability am- 
plitudes. The comp arison with a standard scintillation model 
dBeckert et al.| [2002) indeed suggests good agreement between 
a change of scintillating source size (from 0.07 mas to 0.2 mas) 
and the observed change in variability amplitude (Fig. [3). 

From the model fits of the VLBA data, one can obtain 
an upper limit to the size of the VLBI core and its flux den- 
sity. Thus, one can calculate a lower limit to the brightness 
temperature of that feature. The redshift-corrected lower limit 
brightness temperature (7g ~ (1 + z)S /(v 2 8 2 )) of the core in 
our observations at 5 GHz is > 8.8 x 10 n K. If the brightness 
temperature of the core is equal to this lower limit, assuming 
10 12 K for the inverse-Compton Umit of the brightness temper- 
ature (Kellerm ann & Pau linv-Tothl ll969l) . we find no need for 
relativistic beaming. However, according to earlier VLBA data 
from the archive (year 2004), the redshift-corrected brightness 
temperature of the core at 5 GHz was > llxl0 12 Kat that 
time. Here, at least a Doppler factor of ~ 1 1 is required to re- 
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Fig. 7. Flux density measurements of J 1128+592 by single dish 
and VLBA. Circles represent the single dish flux density mea- 
surements. Diamonds, squares, and triangles show the summed 
model components of VLBA observations at 5 GHz, 8 GHz, 
and 15 GHz respectively. In the insert, all of our 5 GHz sin- 
gle dish observations (circles) are displayed together with our 
5 GHz VLBA data (diamonds), and archival 5 GHz VLBA data 
(crosses). 



duce the brightness temperature to the 10 12 K inverse-Compton 
limit. This means that the minimum Lorentz-factor of the jet is 

Train = <5min/2 > 5.5. 
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Fig. 5. 8 GHz VLBA maps of Jl 128+592 from the six observing epochs. Contours show total intensity, and lines indicate the 
strength and direction of the polarized intensity. Contours are in percent of the peak flux and increase by factors of two. The peak 
flux intensities are 300, 289, 284, 294, 284, and 283mJy/beam. The beam-size is 1.5 mas x 0.7 mas and is shown in the bottom 
left-hand corner of each image. The 1 mas length of the superimposed polarization vectors corresponds to 2.5 mJy/beam. 



To search for possible jet motion on mas-scales, we checked 
for systematic changes in the relative separation of the two VLBI 
components of Jl 128+592. At the highest available resolution, 
at 15 GHz, we did not find any evidence of motion. A linear fit 
to the 15 GHz data formally yields an angular separation rate of 
0.03+0.08 mas/year. This translates into an upper limit to the ap- 
parent speed of (2+6) c. At 15 GHz, we observed a slight change 
in the position angle of the jet feature relative to the core, which 
was assumed to be stationary: in the first three epochs, the po- 
sition angle is ~ -111°, while in the last three it is ~ -107°. 
Compared to the average measurement error in the position an- 
gle of ~ 3°, we regard this change as insignificant. 

We also checked whether any significant motion can be de- 
tected between the 5 GHz archival data from 2004 and our ob- 
servations (2007/2008). Resolving the VLBI structure of the two 
components in the archival data was more difficult because of 
the poorer uv-coverage (Jl 128+592 was a calibrator, usually 
observed in only a few scans.) Therefore, the derived distances 
(and position angles) between the components exhibited a larger 
scatter. A linear regression to the 5 GHz data (weighting of the 
data by the inverse square of the errors) infers a formal speed of 
(-0.02 + 0.02) mas/year, which is consistent with either station- 
arity or marginal inward motion (Fig. [8]). 

Although no significant motion is seen in the jet, we can 
estimate the expansion velocity from the observed change of 
the core size, and thus estimate the viewing angle of the unre- 
solved part of the jet relative to the observer. The core's FWHM 
changed by 0.13 + 0.07 mas during four years, that is at a 
rate of 0.03 + 0.02 mas/year. This translates into an apparent 
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Fig. 8. Relative separation of core and jet component of 
Jl 128+592 as measured with VLBI at 5 GHz. The lines are fits 
from a linear regression analysis. The solid line represents a fit 
with linear weighting of the errors, the dashed line shows a fit 
with quadratic weighting. 

velocity of (2.5 + 1.4)c. The viewing angle is obtained from 
i// — arctan , ■ Using the Doppler factor calculated from 

the brightness temperature of the core component (11), we ob- 
tain for the viewing angle \j/ = 2.3° ± 1°. 
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Fig. 6. 15 GHz VLBA maps of Jl 128+592 from the six observing epochs. Contours show total intensity, lines show the strength, and 
direction of the polarized intensity. Contours are in percent of the peak flux and increase by factors of two. The peak flux intensities 
are 232, 229, 220, 236, 232, and 233 mJy/beam. The beam-size is 1 mas x 0.3 mas and is shown in the bottom left corner of each 
image. The 1 mas length of the superimposed polarization vectors corresponds to 2.5 mJy/beam. 



In summary, we note that J 1 128+592 is at least a mildly (y > 
5.5) relativistic source, which is oriented almost along the line of 
sight. It therefore shows a relatively short jet (seen in projection) 
and a relatively large Doppler-factor of 5 > 11. 

4.1. Annual modulation model and VLBA source structure 

The orientation of the VLBI structure of J 1 128+592 is very sim- 
ilar to the orientation of the scattering ellipse deduced from the 
anisotropic annual modulation model. The average position an- 
gle of the jet feature is 116° at 5 GHz and -108° + 2° at 

15 GHz (note, at 15 GHz, that the position of the jet component 
is located more accurately owing to a smaller beam size). While 
the position angle of the scattering ellipse is -98° ±5°, according 
to the annual modulation model. 

Using the results of our VLBA observations, we fitted the 
variability timescales again, this time keeping the position an- 
gle of the anisotropy fixed at the position angle of the jet fea- 
ture (-116°). The model parameters obtained agree well within 
the errors with the original unconstrained fit (see Sect. |2j. This 
may indicate that the anisotropy seen in the scintillation pattern 
originates mostly in the source structure rather than the scat- 
tering medium. However, additional observations are required 
to clarify unambiguously the roles of the scattering plasma 
and source structure in the anisotropic scattering. According to 
Rickett et al. (2002), a negative "overshoot" in the light-curve 
autocorrelations is indicative of anisotropy caused by the scat- 
tering medium. The study of the scintillation in all Stokes pa- 
rameters (/, U, Q, V) can also shed light on the origin of the 
anisotropy. This analysis however must await future single-dish 



observations of the polarized intensity of J 1128+592, since the 
Urumqi telescope can only perform total intensity observations. 

As mentioned in Sect. [2] annual modulation models of differ- 
ent IDV sources always invoke anisotro pic scattering . The ratio 
of anisotropy is usually above 4 (e.g., Bignall et al. 20 06j re- 
port an axial ratio for PKS1257-326 ~ 12). Walker et al. (2009) 
studied whether an even higher degree of anisotropy can be 
used to describe the observations in the case of J1819+3845 
and PKS 1257-326. The authors found that the anisotropy can 
be so large in these two sources that the scintillation pattern can 
be more effectively described one dimensional. Therefore, they 
concluded it would be more likely that the anisotropy is caused 
by the scattering material rather than the source intrinsic struc- 
ture. 

We used the source model from the VLBI data (created by 
DIFMAP) to investigate how the velocity vector cuts through the 
source structure at the different observing epochs. We used an 
"average" source structure with a core-component flux density 
of 0.26 Jy, a core-component size of 0.2 mas, a jet-component 
flux density of 0.12 Jy, and a jet-component size of 1.1 mas. The 
position angle of the jet component was -116° and its distance 
from the core was 0.5 mas. In 2007.774 (VLBA epoch C), the 
source structure was aligned almost parallel to the velocity vec- 
tor. The direction of the velocity vector at that time of the year 
is -126° (see the ellipse in Fig. fTTb. We cut through the source 
structure in this direction and fitted the resulting curves with a 
Gaussian. The resulting FWHM was 0.17 mas. 

Assuming this to be the value for the major axis of the scin- 
tillating source 6< ma j and using a scintillation length-scale s and 
axial ratio r from the annual modulation model, we can calcu- 
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Fig. 9. The difference between the core and jet EVPA at all three 
observing frequencies measured for the six VLBA epochs. Filled 
black squares represent 5 GHz data, open black squares repre- 
sent 8 GHz data, and gray dots represent 15 GHz data. 



late the screen distance D by adopting the equation fl ma j = 

s ■ (2.54)" 1 • y/r = f ma j ' D - The constant factor of 1/2.54 is 
ca used by a different de finitions of the variability timescale used 
bv lBignall et al.l d2006l) and ourselves. More precisely, since the 
VLBA model fits provide only an upper limit to the source size, 
the derived screen distance represents a lower limit. As a lower 
limit to the distance of the scattering screen, we obtain ~ 37 pc. 

This value is comparable to those derived for fast scintilla- 
tors ( < 10 pc), dDennett-Thorpe & de Bruvnll2003t iBignall et ail 
l2006h . is 3 - 4 times lo wer than the lower lim its obtained previ- 
ously fo r J 1 128+592 in lGabanvi etail d2007bl) . In lGabanvi etafl 
d2007bl) . we used a similar anisotropic annual modulation model, 
but we estimated the source size from the frequency dependence 
of the variability strength at that time. This discrepancy can 
be explained by a not-so-precise estimate because of the lim- 
ited number of simultaneously measured multi-frequency light- 
curves and possible opacity effects. 



5. Results of the VLBA observations - polarized 
intensity 

Both the core and the jet component are polarized at 5 GHz and 
8 GHz. The polarization of the jet component at 15 GHz can be 
detected above the noise level in only three epochs (in 2007.54, 
2007.907, and 2008.025). 

The peak of the total intensity and the polarized feature in the 
core coincides in every epoch and frequency, except in 2007.54 
at 15 GHz. In this data set, the polarized intensity feature is dis- 
placed by ~ 0.15 mas to the east from the total intensity peak. 

The polarized flux densities and the EVPAs of the core and 
the jet are summarized in Tables[3]S and|5] At 15 GHz, we show 
values corresponding to the polarized patches at the largest core 
separation. 

In Fig. |9j we show the difference between the core and jet 
EVPA. Note that this quantity is not affected by a possible cali- 
bration error of the EVPA absolute orientation. (EVPA miscali- 
bration might occur in at least one epoch at 8 GHz, see Sect. [3]). 
It is clear from Fig.|9j that the relative position-angle difference 
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Fig. 10. The rotation measure of the core (filled circles) and the 
jet (open circles) of Jl 128+592 during the VLBA observations. 
The VLBA epochs are denoted by letters (see Table[TJ. 

Table 6. RM of the core (Col. 3) and jet component (Col. 4). 



Epoch 


RMcoie 


>ad/nr] 


tfM jet [rad/m 2 ] 


A 


2007.540 


-501 


± 107 


-65 ± 16 


B 


2007.661 


-90 


±29 


(83 + 159)1 


C 


2007.773 


15i 


:46 


(121 ±46)" 


D 


2007.907 


30 


1 1 


181 ±28 


E 


2008.025 


-48 


±51 


128 ± 127 


F 


2008.137 


-141 


±23 


(86 ±61)" 



" The values in parentheses are based upon only two frequency mea- 
surements. 



between the electric vectors of the two components varies sig- 
nificantly from epoch to epoch (and thus on timescales of ~ 6 
weeks) at all three frequencies. At 5 GHz, the relative EVPA of 
the two features changes by ~ 50°, then it rotates back by 20°. 
At 8 GHz, the trend is similar ~ 60° rotation during the first 
three epochs, then the EVPA difference decreases by more than 
30° to around 100° in the last epoch. In the following, we dis- 
cuss these simultaneous change in the EVPA of two physically 
distinct VLBI components. 

5. 1 . Discussion of the EVPA changes 

In general, changes in the EVPA with frequency are explained by 
the rotation measure of an ionized plasma. When polarized light 
passes through a magnetized plasma, it is subject to Faraday ro- 
tation: the intrinsic polarization angle xo rotates by an a mount 
propo rtional to the square of the observed wavelength dBurnl 
1966). Thus, the observed polarization angle is x — Xo + RM/l 2 . 
The RM depends on both the free electron density n e of the 
plasma along the line of sight and the line-of-sight component 
of the magnetic field Bf RM ~ J n e B\\ds. 

For the jet component, we detected polarization at 15 GHz 
in only three epochs, and we can therefore only obtain a reliable 
RM value for these epochs. We note that plotting the EVPA ver- 
sus lambda 2 , all data points lie along a straight line with relative 
little scatter and are therefore aligned without requiring 180 de- 
gree ambiguity corrections. We also estimated the RM for the 
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Table 3. The results of the 5 GHz VLBA observations. 



Epoch 




Cors 








Jet 








7 C [mjy] 


P c [mJy] 


Xcl° 




Zj [mJy] 


Pj [mJy] 




Urc -x\\l°] 


2007.540 


286 ±5 


3.2 ±0.3 


-24 ± 


4 


89 ±4 


1.2 ±0.1 


83 ±4 


107 


2007.661 


239 ± 17 


3.2 ±0.2 


-30 ± 


4 


134 ± 13 


1.0 ±0.1 


107 ±4 


137 


2007.773 


268 ±8 


4.5 ±0.1 


-49 ± 


4 


93 ±6 


0.9 ±0.1 


96 ±3 


145 


2007.907 


236 ±8 


5.0 ±0.1 


-37 ± 


3 


118 ±6 


0.7 ±0.1 


99 ±5 


136 


2008.025 


229 ±8 


4.2 ±0.2 


-23 ± 


4 


121 ±7 


0.8 ± 0.2 


106 ±4 


129 


2008.137 


243 ±6 


4.3 ±0.2 


-25 ± 


5 


117 ± 9 


0.7 ±0.1 


103 ±5 


128 



Col. 1 lists the observing epochs; Col. 2 lists the integrated flux density of the core component obtained from DIFMAP model-fitting; Col. 3 
lists the polarized flux density of the core component; Col. 4 lists the EVPA of the core; Col. 5 lists integrated flux density of the jet component 
obtained from DIFMAP model-fitting; Col. 6 lists the polarized flux density of the jet component; Col. 7 lists the EVPA of the jet; and Col. 8 lists 
the difference between the EVPA of the core and the jet. 

Table 4. The results of the 8 GHz VLBA observations. 



Epoch 




Core 






Jet 








k [mJy] 


P c [mJy] 


Xcl°] 


Zj [mJy] 


Pj [mJy] 


X\H 


Urc -*jl[°] 


2007.540 


285 ± 10 


2.0 ±0.2 


12+17 


75 ±21 


1.7 ±0.1 


87 ± 17 


75 


2007.661 


264 ± 10 


3.8 ±0.4 


-23 ±21 


81 ±8 


1.2 ±0.2 


96 ±21 


119 


2007.773 


252 ±4 


4.8 ±0.1 


-56 ±4 


93 ± 6 


1.2 ±0.2 


80 ±4 


136 


2007.907 


264 ± 10 


4.1 ±0.2 


-41 ± 3 


84 ± 6 


1.0 ±0.2 


82 ± 15 


123 


2008.025 


256 ±3 


4.7 ±0.3 


-20 ±4 


89 ± 5 


1.1 ±0.1 


83 ±3 


103 


2008.137 


264 ± 10 


4.5 ±0.1 


-9 ±4 


75 ± 10 


1.0 ±0.1 


92 ±3 


101 



The columns are the same as in Table[3] 
Table 5. The results of the 15 GHz VLBA observations. 



Epoch 




Core 






Jet 








I c [mJy] 


Pc [mJy] 


Xd°] 


Zj [mJy] 


Pj [mJy] 


x\n 


Urc -xjlll 


2007.540 


257 ±2 


1.3 ±0.3 


75 ± 10 


39 ±4 


0.5 ±0.3 


96 ± 10 


24 


2007.661 


250 ± 20 


2.7 ±0.1 


-11 ±6 


35 ±4 








2007.773 


238 ± 20 


3.0 ±0.1 


-48 ±5 


26 ±4 








2007.907 


261 ± 11 


2.3 ±0.1 


-42 ± 10 


33 ±3 


0.6 ±0.3 


64 ± 10 


106 


2008.025 


267 ± 10 


3.5 ±0.1 


-4± 10 


36 ±2 


0.5 ±0.1 


113 ± 10 


114 


2008.137 


251 ±3 


2.9 ±0.1 


2±5 


31 ±2 









The columns are the same as in Table[3] 



other epochs with polarization detections at only two frequen- 
cies. Naturally, these values are only approximate estimates of 
the RM (shown in parenthesis in Table |6}. However, it can be 
seen from Fig. [10] that the RM changes smoothly from epoch to 
epoch. 

The RM values obtained are displayed in Fig. [10] The RM 
of the core varied from ~ -500rad/nr to ~ +30rad/m 2 , then 
gradually changed back to a negative value of ~ -141rad/m 2 . 
The RM of the jet also exhibited variations, although less pro- 
nounced. Despite the relatively large measurement errors, it ap- 
pears as if the RM of the core and the jet vary in unison. 

Temporal variability of the RM has been reported for 
only a few AGNs and on timescales of typically 1-3 
years: 3C273 (lAsadaet al.ll2008tlZavala & Tavloril2001l) . 3C279 
(Zavala & Tayl or] |200lL and references therein) and 3C120 
dGomez et alj |200&). The Faraday screen responsible for the 
rotation of the EVPA in these sources is understood to be 
source intrinsic, i.e., either circumnuclear (related to the nar- 
row line region or the torus surrounding the AGN), or internal, 
meaning that it is physically as sociated with the jet (e.g., re- 
lat ed to a sheath aroun d the jet, llnoue et aLll2003l) . According 
to iGomez et alj d2008l) (and references therein), the internal 
Faraday rotation causes depolarization in the source. Thus, when 



the rotation reaches an angle of ~ 45°, the fractional polariza- 
tion decreases by a factor of two. This effect is not seen for 
J 1128+592 in our data, therefore an external Faraday screen is 
more plausible. 

From Fig. [TO] it is obvious that the apparent variation in 
the RM of core and jet exhibit similar trends, despite their pro- 
jected linear distance of > 5 pc. This suggests that the core and 
jet are both covered by the same Faraday screen. The relative 
offsets between the RM of the jet and core, and the slightly 
more pronounced variability in the RM of the core component 
is indicative of a second Faraday screen, which covers only the 
core region but not the jet. We tried to decompose the RM into 
two Faraday screens, one that is common to both core and jet, 
and causes a similar variability trend throughout our six VLBA 
epochs (ARM ~ 150rad/m" on a timescale of 5 to 6 months), 
and a second component, which only affects the core-region of 
the VLBI structure. This latter Faraday screen causes an RM 
variation in the core-region of ARM < 300rad/m 2 on shorter 
timescales of 1.5 months. At this point, it is obvious that this 
decomposition is affected by the relatively large measurement 
uncertainties and should be confirmed by future more finely (in 
time and frequency) sampled measurements. 



10 



K. E. Gabanyi et al.: Very Long Baseline Array observations of the intraday variable source Jl 128+592 




Fig. 11. The variation in the rotation measure. The ellipse (solid 
black line) shows the apparent velocity of the screen relative to 
the Earth in the two orthogonal coordinates of right ascension 
and declination according to the fitted anisotropic annual modu- 
lation model. The dashed lines (blue) represent the velocity vec- 
tors at the times of the 6 VLBA observations (labeled A to F, see 
Table[TJ. The vertical lines in the plane orthogonal to the orbital 
motion indicate the absolute value of the RM at the times of the 
VLBA observations. 



We may hypothesize that the screen affecting both compo- 
nents identically could be the same local screen that is respon- 
sible for the IDV and the annual modulation. In this case, the 
observed systematic variation in the RM could be caused by: 
(i) variations in the parallel component of the interstellar mag- 
netic field (ii) variation in the electron density n e in the local 
screen, or (iii) a variation in the path length, i.e., the thickness 
of the screen with respect to its projection on the ecliptic (or- 
bital plane), or of course, a combination of these three effects. If 
we attribute most of the observed systematic RM variability to a 
screen in the nearby galactic ISM, changes in the magnetic field 
and/or the electron density over spatial scales of shorter than 2 
AU (or timescales of a few weeks) would indicate the presence 
of MHD-turbulence in the ISM. On the other hand, if the mag- 
netic field is homogeneous on AU-scales and frozen in, it will 
not vary strongly (with neither time nor on AU-size scales). The 
projection angle of the magnetic field (B\\ ~ Bq sin(</>)), could 
however in principle vary if the screen is very close, because the 
Earth orbits the Sun. For a presumed screen distance of ~ 37 pc, 
one expects variations in the projection angle of no more than 
~ arcsin(2AU/37pc) ~ 10~ 5 °, which is far too small to cause 
any significant variation in the RM. It is also difficult to ex- 
plain a fractional RM variation of A RM ~ 150 rad/m 2 solely by 
changes in the path length. This points towards a combination 
of time-variable path length, electron density, and perhaps even 
magnetic field, in such a way that it would explain the observed 
systematic RM changes, as illustrated by the 3D-diagram in Fig. 
[TT1 We note that this time dependence would be the natural con- 
sequence of a refractive lens moving through the line of sight. 
These lenses are used to explain t he so-called "extreme scatter- 
ing events" in some IDV s ources dClegg et alJll998tlLazio etail 
2000; Senkbeil et al]|2008l) and cause systematic in tensity vari- 
ations on timescales of months. IClegg et all dl998l) discuss the 
geometrical effects of such a lens and mention the possibility 



of RM variations, which, in their case, were not however de- 
tected. Since the frequency-dependent "bending" of the rays by 
the plasma lens leads to caustics and differential amplification of 
the (polarized) substructure of the source, it is possible that such 
a lens also introduces apparent RM variations. The discussion of 
this possibility, however, is beyond the scope of this paper and 
deserves a future, more thorough elaboration. The variation in 
the line-of-sight "thickness" of the screen (or the related emis- 
sion measure, which is ~ J n e ds) should also produce variation 
in the scattering measure (the line-of-sight integrated rms elec- 
tron density) and the related variability index. However, the lim- 
ited accuracy of the variability-index measurements illustrates 
these variation s only very mar ginally (see solid circles in Fig. 3). 
According to lSpangrerld2009l) . the RM of ISM clouds located in- 
side the Local Bubble is typically of the order of 1 - 10rad/m 2 . 
This is much lower than the observed RM variations towards 
J 1 128+592, or even the part of RM variations, which is the same 
in both jet and core according to the decomposition analysis. 

At this point, the signal-to-noise ratio of our data (i.e., RM, 
IDV timescale, and IDV amplitude) is clearly insufficient to dis- 
entangle the RM variations caused locally (interstellar weather, 
line-of-sight effects) and in the source itself. Higher polarization 
sensitivity, a more accurate determination of the RM of VLBI 
core and jet, and a longer time coverage (several orbital periods) 
will be necessary to study this effect further. 

In an alternative and perhaps more conservative approach, 
one can argue that we do not see an intrinsically variable RM, 
nor does the intrinsic EVPA of the VLBI components change 
with time: we are observing instead at different times blends of 
different subcomponents, which cannot be resolved in our VLBI 
observations, even at 15 GHz. This blending of polarized com- 
ponents were reported in the case of other IDV sources, e.g., 
PKS0405-385 dRickett et al.ll2002l) . Thus, the EVPA changes in 
J 1 128+592 might also be explained by the emergence of a new 
polarized feature. As the new polarized subcomponent gradu- 
ally separates from the core, the core EVPA changes, since it 
is dominated more and more by the new, strongly polarized 
feature. As this feature expands and thus fades, the measured 
EVPA becomes dominated again by the original core polariza- 
tion, and thus it rotates back towards its original value. This 
scenario would naturally explain the observed shift between 
the total intensity and polarized intensity features in the first 
epoch at 15 GHz. (Similar offsets are often reported in quasars 
and blazars, and explained as polarized subcomponents e.g. 
iBach et alJbOOfjt iGabuzda & G6mezll200Tl) . However, this sce- 
nario would not explain the EVPA variations (and the similar 
trend in the core and jet RM variations) of the jet component. 
Moreover, we also observed a variation in the RM of the core and 
jet without significant time shift. This would be difficult to un- 
derstand, if the RM of the core were affected by blending effects, 
e.g., because of a newly emerging polarized jet component. 

A solution to this problem and a viable interpretation of the 
correlated variability in the EVPAs and the RM in core and 
jet, must most likely include two effects: (i) a Faraday screen, 
which covers the entire VLBI structure and leads to the observed 
similar variability trend in the RM of jet and core, and (ii) an- 
other process, which only affects the frequency alignment of the 
EVPAs in the VLBI core region, and for which only marginal ev- 
idence is seen in our data. At present, it is unclear whether this 
second process is related to blending effects in a time-variable 
polarized substructure of the core or to another Faraday screen 
(in the quasar), which affects only the core but not the jet (cir- 
cumnuclear screen, e.g., a torus). Future polarization sensitive 
VLBI monitoring of higher sensitivity, denser frequency cover- 
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age, and higher angular resolu tion - i.e., as provided by the fu- 
ture space VLBI (VSOP2, e.g. Tsuboi 2008) - should help us to 
distinguish between source-extrinsic and source-intrinsic depo- 
larization effects. 



6. Summary 

We have presented new results derived from our 3 frequency 
(5/8/15 GHz) polarization-sensitive VLBA observations of the 
IDV source Jl 128+592. An aim of these observations was to 
discover whether source structural changes are responsible for 
the decrease in the strength of IDV observed in this source dur- 
ing 2007-2008. 

Our data uncovered an east-west oriented jet-core structure 
of the source. For the entire six epochs of VLBA data (span- 
ning 0.6 yrs), we were unable to detect any significant change 
in the total intensity VLBI structure of the source. Compared 
with archival VLBA data, we found that the total flux density 
of the core changed significantly. This agrees with our single- 
dish monitoring result, which included a flare-like brightening 
in Jl 128+592 during early 2006. Compared to archival data, the 
size of the core at 5 GHz also increased significantly, which may 
explain the quenching of the intraday variability. Although we 
did not detect any structural changes in the source during our 
observations, the core size was found to be excessively large at 
15 GHz. This may be indicative of the emergence of a new jet 
component that is still blended with the core at the time of our 
observation. 

We fitted an annual modulation model to the IDV timescales 
derived from our densely time-sampled IDV flux-density moni- 
toring data. One important parameter of the annual modulation 
model is the scattering lengthscale, which is determined by the 
scintillating source size and distance to scattering screen. From 
the IDV monitoring and the size limits of the VLBI core, we ob- 
tained a lower limit to the distance of the scattering screen of 
~ 37 pc. 

The source structure detected by our VLBA observations has 
an orientation that roughly agrees with the orientation of the 
anisotropic scattering obtained from the fit of annual modula- 
tion model. This might indicate that the source-intrinsic struc- 
ture is responsible for the anisotro pic scattering (unl ike in other 
well-known IDV sources, e.g., Walk er et al.l 120091 and refer- 
ences therein). 

We found significant variations in the orientation of the po- 
larization vector in the source (core and jet) with time. At all 
three observing frequencies, the EVPA and the deduced RM of 
the core varied in unison on timescales as short as 6 weeks. 
The variation in the RM of the jet is less pronounced than in 
the core and affected by larger measurement errors, but seems 
to follow the same systematic trend as the RM variability of 
the core. The relatively fast variation of the RM and the simi- 
lar trend in RM variations of the core and jet might indicate that 
at least a certain fraction of the RM variability is caused by a 
nearby ISM screen, which can also be responsible for the IDV 
and the annual modulation of the variability timescales. At this 
point, other source-intrinsic interpretations of the RM variability 
(e.g., internal Faraday rotation, different RM in different regions 
of the source) cannot be excluded. More data with superior time 
and frequency coverage will be necessary to distinguish between 
these possibilities. 
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